Tunable Quantum Fluctuation-Controlled Coherent Spin Dynamics 
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Temporal evolution of a macroscopic condensate of ultra-cold atoms is usually driven by mean 
field potentials, either due to scattering between atoms or due to coupling to external fields; and 
coherent quantum dynamics of this type have been observed in various cold atom experiments. 
In this article, we report results of studies of a class of quantum spin dynamics which are purely 
driven by zero point quantum fluctuations of spin collective coordinates. Unlike the usual mean-field 
coherent dynamics, quantum fluctuation-controlled spin dynamics or QFCSD studied here are very 
sensitive to variation of quantum fluctuations and the corresponding driving potentials induced 
by zero point motions can be tuned by four to flve orders of magnitude using optical lattices. 
These dynamics have unique dependence on optical lattice potential depths and quadratic Zeeman 
fields. We also find that thermal fluctuations generally can further enhance the induced potentials 
although the enhancement in deep optical lattices is much less substantial than in traps or shallow 
lattices. QFCSD can be potentially used to calibrate quantum fiuctuations and investigate correlated 
fluctuations and various universal scaling properties near quantum critical points. 



I. INTRODUCTION 

When particles such as atoms interact with each other 
at low temperatures, very often they exhibit remarkably 
distinct cooperative behaviors as a result of symmetry 
breaking. One of the most fascinating and distinct con- 
sequences is the possibility of observing quantum dynam- 
ics at a macroscopic level [H, In Bose-Einstein con- 
densates of ultra-cold alkali atoms, macroscopic quan- 
tum phenomena related to coherent matter waves [j, 0], 
AC Josephson effects [1, 0], and vortex lattices Q have 
all been observed. Studies of these phenomena in ultra 
cold matter should eventually lead to applications such 
as cold- atom-based precise measurements. 

Spin correlated macroscopic quantum dynamics have 
also been a focus of many cold atom experiments car- 
ried out recently. Spin ordering and spin-relaxation col- 
lisions were first investigated in condensates of sodium 
atoms by Inouye et al. and Miesner et al. from Ketterle's 
group at MItH, i]. Coherent spin dynamics driven by 
various mean field interactions or external fields were 
later demonstrated in condensates of hyperfine spin-two 
rubidium atoms (l0|. and hyperfine spin-one rubidium 
atoms [ni Ell ■ Ordering in spinor gases is usually induced 
by hyperfine spin dependent two-body scattering [isl [l^ . 
Coherent spin dynamics observed in experiments are re- 
lated to the coherent quantum dynamics explored in solid 
state superconductors, and earlier experiments on ultra 
cold gases of atoms [a, @. They are explicit manifes- 
tations of fascinating macroscopic quantum states and 
can be potentially applied towards constructing high pre- 
cision interferometers. Remarkably, coherent dynamics 
also provide a unique direct measure of interaction ener- 
gies or scattering lengths as emphasized before 5, 10, 11]. 

Quantum fluctuation-controlled spin dynamics or 
QFCSD we are going to study in this article on the 
other hand are a direct measure of quantum fluctua- 
tions; they can be potentially used to calibrate quantum 
fluctuations and investigate correlated fluctuations near 



quantum critical points or universal scaling properties. 
Furthermore, QFCSD of cold atoms can be designed to 
simulate many other quantum-fluctuation-induced phe- 
nomena such as Coleman- Weinberg mechanism of spon- 
taneous symmetry breaking 15 1, and order due to dis- 
order in antiferromagnets 16r[l3]- The system we are 
examining to understand QFCSD is a condensate of ru- 
bidium atoms {^'^ Rb) in hyperfine spin-two {F = 2) 
states. The two-body scattering lengths between rubid- 
ium atoms have been estimated using both photoassocia- 
tion data J3, 24, 25, 26] , and elastic scattering data near 
Feshbach resonances [23, [28|. Most recently, Rabi oscil- 
lations between different two-spin states have also been 
used to measure spin dependent interactions [29j . 

Quantum fluctuations of spin collective coordinates 
at least have three different effects on spin dynamics. 
Firstly, quantum fluctuations of wave lengths of conden- 
sate sizes usually result in finite-size quantum symmetry 
restoring. For a spinor condensate of a few million ^'^ Rh 
(F=2) or "^^Na (F=l) atoms, restoring of spin rotational 
symmetry or quantum diffusion of spin orientation typ- 
ically occurs at a time scale of a few tens of seconds. 
So they are not relevant to 10 — 100ms spin dynamics 
studied in large condensates. The second well-known ef- 
fect is to renormalize semi-classical spin dynamics, such 
as spin wave velocities. In high dimensions, dominating 
contributions are usually from short wave length fluc- 
tuations (in optical lattices the shortest wave length is 
set by lattice constants). These spin fluctuations might 
also result in quantum phase transitions between spin or- 
dered and disordered states. For most dilute cold atom 
gases, the renormalization of spin dynamics is perturba- 
tive and negligible (except around critical points). And 
in the limit that interests us, the relative amplitude of 
spin fluctuations is small. The third effect is to in- 
duce dynamics when conventional mean field dynamics 
are completely frozen out because of degeneracies in a 
submanifold. This is the limit of QFCSD that we are go- 
ing to study. The dominating contributions in our case 
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are mainly from quantum fluctuations of an intermediate 
wave length which is comparable to the de Broglie wave 
length of an atom traveling with spin wave velocities. 

However it is quite challenging to probe QFCSD in 
traps without optical lattices, if not impossible. For di- 
lute gases in the absence of optical lattices, the effec- 
tive driving potential induced by quantum fluctuations 
is about lO^^pk per particle (see section III for more dis- 
cussions) because the relative amplitude of fluctuations 
is very small. The corresponding dynamics driven by 
such a small potential are only visible at a time scale of 
a few thousand seconds, too slow to be observed in cur- 
rent cold atom experiments. In addition, a tiny external 
magnetic field of strength ImG can result in a quadratic 
Zeeman coupling of order of 6 x 10 ^'^pk for rubidium 
atoms which is three orders of magnitude larger than the 
induced potential in dilute gases. In most experiments 
because of noises in lasers, the effective quadratic Zee- 
man coupling can be controlled only up to a uncertainty 
that is equivalent to a magnetic field of an order of ImG. 
This further complicates future experimental studies of 
QFCSD. To resolve these difficulties, we propose to en- 
hance the effect of QFCSD using optical lattices. To vary 
the amplitude of quantum fluctuations and optimize the 
effect of fluctuations, we study QFCSD in optical lat- 
tices where the optical potential depth y is a convenient 
tunable parameter [H [li, [13, [HI, . 

The rest of the article is organized as follows. In sec- 
tion II, we introduce a lattice Hamiltonian to study dy- 
namics of Rb atoms in optical lattices and discuss the 
range of parameters we have used to investigate this phe- 
nomenon. In section III, we present our main numerical 
results on quantum-fluctuation induced potentials, fre- 
quencies of coherent dynamics, and how potentials and 
frequencies depend on optical lattice potential depth. We 
also study the dynamical stabilities of coherent oscilla- 
tions when a quadratic Zeeman coupling is present. In 
section IV, we discuss effects of thermal fluctuations and 
analyze the potential-depth dependence of thermal en- 
hancement of induced driving potential. In section V, 
we further investigate effects of spin exchange losses and 
propose how to observe quantum-fluctuation-controlled 
spin dynamics within a relatively short life time of = 2 
rubidium atoms. In section VI, we conclude our studies 
of QFCSD. 
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Here k is the lattice site index and < kl > are the nearest 
neighbor sites, fi is the chemical potential and tL is the 
one-particle hopping amplitude, qb is the quadratic Zee- 
man coupling constant, ql, and cl are three interac- 
tion constant which have been calculated (see subsection 
B for more discussions). The single band Hamiltonian 
is valid when all interaction constants above are much 
smaller than the energy spacing between centers of two 
lowest bands. 

We have employed the traceless symmetric matrix op- 
erator ?/;^ that was introduced previously for the studies 
of hyperfine spin two rubidium atoms '3Cr|; components 
^1/3' = 2/j -2 are linear superpositions of five spin-2 
creation operators, V^mpj — 0, ±1, ±2 as given above. 
It is advantageous to use this tensor representation if 
one is interested in rotational symmetries of condensate 
wavefunctions, or construction of rotationally invariant 
operators. We use it to analyze collective spin modes 
that correspond to small rotations around various axes. 
The tensor operator ■i/)'!' is defined in terms of the usual 
creation operators i/ij^ for five F = 2 states. 



mF=0,±l,±2 



C,,(±l) = iCy(T2) 



72' 



V2' 



CzM 



V3 



(2) 



where Capimp) is symmetric with respect to a/3; all other 
coefficients are zero. The number operator p, the dimer 
or singlet pair creation operator I?^ , the total spin opera- 
tor Fa are defined as p = 1 / 2trip^ ip , = 1 / ^/40trilj^ 
Fa — — ieQ/37V'^^'0r;7- And the quadratic Zeeman opera- 
tor Qzz is defined as Qzz = tri/j^Qtp, and Qap = SazSpz- 



B. Range of interaction parameters for the lattice 
Hamiltonian 



II. MODEL 

A. Microscopic Hamiltonian for spin-two rubidium 
atoms in optical lattices 

In optical lattices, we use the following Hamiltonian 
that was introduced previously [sol [s^ . 

k 



Spin correlations between hyperfine spin-two rubidium 
atoms are determined by three two-body s-wave scat- 
tering lengths Of, F = 0,2,4. In optical lattices, local 
spin-dependent interactions contain two contributions as 
shown in Eq.(IT]); one is, feL-^|/2, the energy of having 
total hyperfine spin J^k at site k, and the other is the 
energy of creating spin singlet pairs (dimers), bcL'Dl'Dk 
where Vk is the dimer creation operator [s^. The usual 
contact interaction at site k is of the form aL{p1 — pk)/2, 
where pk is the number of atoms. Three effective cou- 
pling constants aL,bL,CL which characterize various in- 
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V[E^] V[E^] V[Ej^] 



FIG. 1: (Color online) Coupling parameters ul, bL, cl and 
hoping integral t_t (all in units of nk) as a function of optical 
potential depth V (in units of recoil energy Er) in 3D (black 
circles) and 2D (red squares) optical lattices. Er = 3.18^7? 2: 
for A — 850nm lasers. 



FIG. 2: (Color online) Condensation fraction versus optical 
lattice potential depth V (in units of i5_R).The inset shows the 
condensate fraction at different temperatures when V is equal 
to 6Er. 



teractions are functions of two-body scattering lengths 
ap, F = 0,2,4: and on-site orbitals V'o(r), 



aL{hL,CL) = a{b,c) 



dv[%{v)M^)f- (3) 



Here a = (4a2 -|- 3a4)/7, h — (04 — a2)/7 and c = (7ao — 
10a2 -I- 3a4)/35 are three effective scattering lengths; 
is the localized Wannier function obtained by solving the 
Schrodinger equation for an atom in periodical poten- 
tials. aL,bL,CL can then be calculated usin g th e esti- 
mates of scattering lengths obtained in Ref. [2g. The 
range of these parameters for rubidium atoms in optical 
lattices is plotted in Fig[TJ 

The range of lattice potential depth is chosen to be 
from zero up to ten recoil energy or lOEp;. {Er = 
h? /2mX^, where A is the wavelength of lasers.) Below 
iEji, quantum fluctuations turn out to be too weak to in- 
duce substantial potential and dynamics. For this reason, 
we only show numerical results for V larger than 2Eft but 
less than IGEr- Within this range, we find that quan- 
tum depletion is usually small, less than twenty percent 
and our lowest order calculation should suffice. Although 
quantum fluctuations can be further enhanced above 
IOEr and a superfluid-Mott phase transition should take 
place at 13-Eb jisl. [l9|. close to a critical region we how- 
ever expect a perturbative calculation like the one carried 
out in this article becomes invalid. The range of atom 
number density here, or the number of atoms per lattice 
site M is from zero to three. Most of data are shown for 
typical values M = 1.0, 2.0, 3.0. 

For F = 2 rubidium atoms, the quadratic Zeeman cou- 
pling qb is related to a uniform external magnetic field 
B via qs = 3{fiBB /2)'^ / A, A = 6.8GHz. Here hb is the 
Bohr magneton and A is the hyperfine splitting. For the 



purpose of studying QFCSD, we set the range of B to be 
TOmG > B > IniG where effects of quantum fluctuations 
are most visible; beyond 70mG the dynamics are mainly 
driven by the mean field quadratic coupling qs (see sec- 
tion VI for more discussions). The corresponding range 
for the quadratic Zeeman coupling is then from lO^^pk 
to lOOpfc. 



C. Mean field ground states: quantum spin 
nematics 

The coherent spin dynamics of a uniform condensate 
can be described by the evolution of a condensate wave- 
function Xj i-G. the expectation value of matrix operator 
■0'!'. The corresponding equation of x is 



1 dx 

2 dt 

dx* 



dx* 



OH,,, dH, 



Q 



dH, 



if. 



dx* ' 



dx* dx* 
= ^xTr(x*x) + YX*Tr(xx) 

+ Y^x,[r,x]]-iztL + ^^/2)x, 

= QbQx- 



(4) 



Here Hsc is the semiclassical Hamiltonian (of matrix x, 
x) obtained previously [30l [3]] . [3^ . Hq is the quadratic 
Zeeman coupling term with qB being the coupling 
strength. H^c and i/g (per lattice site) are given by: 



Hsc = 'Y^'^ix*x)T'^ix*x) + yTr(x*x*)Tr(xx) 
+ ^Tr[r,x]'-(^iL+M/2)Tr(rx); 



Hq = teTr(x*gx). 



(5) 
(6) 
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z is the coodination number of optical lattices. For 
a field along the z-direction, matrix Q is defined as 
Qq/3 = ^az^fiz- And Hqj is the quantum- fluctuation- 
induced Hamiltonian discussed below. 

For rubidium atoms, scattering lengths estimated in 
Ref. [23. [2R HI] lead to = — IOcl while results in 
Ref . I27I, l28'| yield = —2.8cl. Both calculations show 
that interaction parameters satisfy cl < and 46^ > cl. 
As pointed out before (30|. in this parameter region with- 
out quadratic Zeeman coupling, ground states are spin 
nematics characterized by real and symmetric tensor 
wavefunctions x (up to an overall phase). Any conden- 
sate that is initially prepared in this submanifold has 
no mean field dynamics because the potential gradient 
dHsc/dx* vanishes. 

To highlight the structure of degenerate ground states, 
we consider an arbitrary condensate amplitude x in the 
spin nematic submanifold. It can be parameterized using 
an SO {3) rotation, TZ and a.U{l) phase shift (p, and a real 
diagonal traceless matrix x; i.e., 

n{x^,Xy,x,) = cxp(r"x, + T^Xy + t'x,) (7) 

where M is the number density or average number of 
atoms per lattice site. TZ is an 50(3) rotation matrix 
defined by three spin angles Xa, a = x,y,z, and anti- 
symmetric generators Tg, 



Pi 



xiC): C £ [0, 27r], 

are normalized real diagonal traceless matrices that form 
a family of solutions specified by a single parameter f (33|: 



U(l) X SO(3) 




^=0 jr/6 n/3 



FIG. 3: (Color online) An artistic view of the five-dimension 
manifold of spin nematics. Three 50(3) rotation degrees and 
I7(l)-phase degree are represented by a cross section at a given 
spin deformation angle ^. We have also plotted the nematic 
wavefunctions i/'s(6', 0) — Xa/snan/s as a function 9, (j) in 
a spherical coordinate, n^, is the ath component of a unit 
vector n{9, 0) and colors indicate phases of wavefunctions. 



0„/3 = 2Msin(2e + ^c.)^a/3. 



(11) 



Xc 



sin(^ - Ca) 
V3 ■ 



(8) 



and = 7r/6, ^j, = 57r/6 and — 37r/2. Nematics with 
different ^ exhibit different spin configurations. Follow- 
ing the definition of tensor operator in Eq.([2|), one can 
easily show that these solutions represent condensates 
of spin-two atoms specified by five-component wavefunc- 
tions V''^ = (V'2,i/'i,V'o, V'-i, V'-2) and 



' M[ — 0, C0S4, 0, — -j^). 



V2 



V2 



(9) 



Therefore, the nematic submanifold is effectively a five- 
dimension space that is characterized by five collective 
coordinates X^^ v — three spin rotational 

angles Xa (a = x,y,z), one spin deformation ^-angle 
X^ (or ^) specifying spin configurations and one phase 
angle Xp{= (f>). Furthermore, it is easy to verify that 
up to an SO (3) rotation and a phase factor, states at 
^ and at ^ -|- tt/B are equivalent and the fundamental 
period of this characterization is tt/B. When the system 
is rotationally invariant (i.e. no external fields), we also 
find that E{^) = + tt/B) and E{^) = E{-^). E{^) is 
the energy of a state defined by ^. 



So in a more conventional representation, these states 
labeled by value <^ correspond to condensates where all 
atoms occupy a particular spin-2 state. 



D. Effective Hamiltonian for QFCSD 

To investigate the kinetic energy and quantum 
fluctuation-induced potential energy for dynamics of a 
condensate initially prepared in this submanifold, we ex- 
1^ >= cos^|2, > H 7^(|2' > ~2 >). (10) pand tensor i/jl^g about a reference condensate wavefunc- 



sin^ 



V2 

It is worth pointing out that spin nematics here are 
time-invariant; the expectation value of the hyperfine 
spin operator J^^. in these states is zero. However, all 
nematics have the following nonzero quadrupole spin or- 
der (up to an SO (3) rotation). 



a/3 

tion x(C) in terms of five collective coordinates X^ in- 
troduced above and their conjugate operators Pi,. We 
furthermore separate the macroscopic dynamics of con- 
densates (q = 0-mode) from the microscopic zero point 
quantum fiuctuations (q 7^ mode) . And we restrict our- 
selves to the dynamics of a condensate in a linear regime. 
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In the appendix, using the decomposition introduced in 
Eg. ljAip we expand the Hamiltonian in Eq.([T|) in terms 
of collective coordinates X,^, Pi, and obtain an effec- 
tive Hamiltonian for submanifold dynamics. Up to the 
quadratic order, the Hamiltonian contains two sectors, 
one involving operators of q = 0-mode and the other one 
only involving operators of q 7^ 0. The Hamiltonian for 
the q = sector generates the kinetic energy needed for 
the dynamics along five orthogonal directions. The cor- 
responding effective masses can be expressed in terms of 
scattering lengths ao,2,4, and quadratic Zeeman coupling 
qs- In addition, the expansion of the quadratic Zeeman 
term Hq for q = mode also generates a mean field 
potential Vq. This potential Vq as illustrated below al- 
ways favors a biaxial nematic with ^ = 7r/2. The biaxial 
nematic has dihedral-four {Dihi) symmetries with easy 
axes in the a:^- plane [33|. On the other hand, the sector 
of Hamiltonian for q ^ collective modes contains zero 
point energies of those modes. These energies in gen- 
eral depend on spin configurations or the values of ^. So 
quantum fluctuations of collective coordinates effectively 
induce a ^-dependent potential Vqf in the submanifold; 
this potential alone selects out a unique ground state as 
recently pointed out by Song et al [34I and Turner et 
al [3^. For rubidium atoms with a positive b^, previ- 
ous calculations show that the ground state is a uniaxial 
nematic in the absence of quadratic Zeeman coupling. 
Turner et al also pointed out that thermal fluctuations 
further enhance the amplitude of induced potentials and 
this order-from-disorder phenomenon is robust against 
finite temperatures [33]. 

We now study the dynamical consequences of both 
quadratic Zeeman coupling Hq and zero point quantum 
fluctuations Hqf. For simplicity, we are mainly focused 
on dynamics around 

a) a uniaxial nematic at ^ — Q ov a. condensate with 
rubidium atoms occupying hyperfine spin state |2,0 >; 

b) a biaxial nematic at ^ = tt/2 or a condensate with 
rubidium atoms occupying hyperfine spin state (|2,2 > 
+ |2,-2>)/y2. 

The resultant Hamiltonian for oscillations around a 
state ^ = or 7r/2) can be cast in the following form, 



2NTm„ 
AMqB 



Vq{£„X,) 



cos2(^ + X^) + AMqB J2 «-^a- (12) 



In Eq.(|T2|). Nt is the number of lattice sites and Al is 
the average number of atoms per site. The masses for 
five directions are calculated to be 



1 



2M 



flic 

ul + cl' qsK^-^McL 



(13) 







^2 

^xx(yy) 

\z ~ 'ydzz^ is ^ function 



i^x,y,z,(, are ^-dependent; Kx{y) 
and is independent of ^; = 4(x 
of ^ and is —4/3 when ^ = 0, and 4/3 when ^ — 7r/2; 
Xzz = dxzzldi. GaaiOi a = x,y,z IS & function of ^, 
Gxx,yv{0 = sin'(CT f ), G.. = sin^ ^ E„,^{tX^,qB) 
is the energy of a mode-z^ = x, y, z, ^,p) collective ex- 
citation with crystal momentum q, and is a function of 
parameters ^, X^, and quadratic Zeeman coupling qB- In 
the appendix, we show the general form of this energy ex- 
plicitly for the case of nonzero qB- Only spin modes with 
ly = x,y,z,^ contribute to the ^-dependence of poten- 
tial Vqf. Fluctuations of phase modes are independent 
of parameter ^ or Zeeman coupling qB and are irrelevant 
for discussions of spin dynamics as a result of spin-phase 
separation. Note that Vqf is a function of ^, qB and X^; 
and Vq is a function of ^, qB and X,y, v = x,y, z, ^. 



III. MAIN RESULTS FOR QFCSD 
A. Potentials induced by quantum fluctuations 

We first consider a situation where the quadratic Zee- 
man coupling is absent and the potential Vq vanishes. 
As argued before, generally speaking Vqf{^, 0) = Vg/(^ -I- 
7r/3,0) and Vq/(^,0) = Vqf{-^,0). The explicit form of 
Vqf calculated here indeed is consistent with this gen- 
eral requirement. We also find that the main contribu- 
tion to Vqf is from fluctuations of wavelength l/(mBNVa) 
that is much longer than the lattice distance d^. More 
speciflcally, following discussions in the appendix, the ^- 
dependent induced potential can be written as a sum of 
polynomials of spin-wave velocities, i.e. J^a ^q"*"^ ^^^r d- 
dimension optical lattices. Spin-wave velocities are func- 
tions of CL,bL and II- We numerically integrate over 
all wavelengths and obtain the ^-dependence of Vqf as 
shown in the inset of Fig[4l 

Consequently, we find that in d-dimension lattices 
(d = 2,3), the barrier height Bqf which is defined 
as Vqf(^,Q) — Vq/(0,0), the energy difference between 
^ = 7r/6 and ^ = satisfies the following simple scaling 
function, 



B. 



qf 



\Mcl\' 

t? 



■ ,bL M\cl\ 
-9d( — , — 

CL tL 



(14) 



Here gd{x, y) is a dimensionless function that can be stud- 
ied numerically. This scaling function is either insensitive 
to the variation of y as for d = 3 [32] or independent of 
y as for d — 2. 

We also study the barrier height Bqf as a function of 
V, the potential depth of optical lattices. We find that in 
the absence of lattice potentials or in traps, for a density 
that is equivalent to one particle (M = 1.0) per lattice 
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0.01 



0.00 r 



10.0 



FIG. 4; (Color online) Barrier height Bqf (in units of pK) as a 
function of optical potential depth V (in units of recoil energy 
En). Inset is for the ^-dependence of Vqf, a potential induced 
by quantum fluctuations when V — QEr and M = 2.0. 



site the barrier height is of order of 10~^ to 10~^pfc and 
is negligible in experiments. When the optical lattice 
potential depth V is varied, the barrier height typically 
increases by four or five orders of magnitude. Particu- 
larly, as V increases from 2Eii to lOiJ^, mean field in- 
teraction energies aL,bL,CL vary by less than a factor 
of three; however, for the same range of V , the barrier 
height Bqf varies from 10~^pk to a few pk. The energy 
shift between the uniaxial state at ^ = and biaxial state 
at ^ = 7r/6 is analogous to the Lamb shift observed in 
atoms [3^ ■ In Figdl we show the barrier height Bqj of 
induced potential verus lattice potential depth V. 

In Fig. [5l we further plot the effective potential 14/ / = 
VQ+Vqf as a function of optical potential depth V. When 
V is less than 14, the biaxial nematic with ^ = 7t/2 is 
stable and as V > I4, it becomes locally unstable. For 
^'^Rb atoms with one atom per lattice site (M = 1.0), 14 
is about 5.5Eji {En is the recoil energy of optical lattices) 
when the quadratic Zeeman coupling is lOpfc. Almost 
opposite behaviors are found for uniaxial nematics at ^ = 
0. Values of ^ at which global potential minima are found 
are shown as a function of 14// in Fig. [5Jd. 



B. Oscillations induced by quantum fluctuations 

To understand the dynamical consequences of Vqf , we 
consider coherent dynamics around the uniaxial nematic 
state at ^ = that is selected out by the fluctuation- 
induced potential Vqf or the biaxial nematic state at 
^ = 7r/2 that is favored by the quadratic Zeeman po- 
tential Vq. Particularly we are interested in oscillations 
along the X^-direction around state ^ = or 7r/2. These 
motions correspond to the following time evolution of 
condensates, 




-0.01 



^V=9Ej^ 
<^^V=11E„ 




FIG. 5: (Color online) (a) Effective potential 14// = Vq+Vqf 
(in units of nk) as a function of ^ for various optical potential 
depth V (in units of recoil energy Er) of 3D lattices. The 
quadratic Zeeman coupling strength is set to be lQpk{31mG) 
and number of atoms per lattice site M is equal to one. (b) 
Values of ^ at which global potential minima in (a) are found 
are plotted as a function of optical potential depth V. 




FIG. 6: (Color online) (a),(b), frequencies fo — no/(27r) for 
population oscillations around |2,0 > state are shown as a 
function of potential depth V in (a) 2D and (b) 3D opti- 
cal lattices for different atom number density M; here the 
quadratic Zeeman coupling is absent. 
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\t (cos^-sinCXoF(t))|2,0 > + 
-^(sin^ + cos^XoF(t))(|2,2 > +|2,-2 >); 



Fit) 



cos —t + i—T-r sm —t 
2 4Af 2 



(15) 



where is the initial deviation from state ^ = or 
7r/2 and is the effective mass given in Ea. (IT3)) . By 
solving the equation of motion in Eq. (|12p including the 
quadratic Zeeman effects, we derive a general expression 
for oscillation frequencies ri(= 27r/). The oscillation fre- 
quency il{qB,V) in general is a function of the quadratic 
Zeeman coupling qg, optical lattice potential depth V, 
and the number density M. For population oscillations 
around ^ = or 7r/2, we obtain the following frequencies 



= 2 



8Mqi 



(16) 



V^j is the 



here cr = ^1 for ^ = and 7r/2 respectively, 
curvature of potentials at ^ = or 7r/2. 

Note that both the real and imaginary part of F{t) 
oscillate as a function of time; thus unless -jj^ is exactly 
equal to one which only occurs when atoms are noninter- 
acting or the quadratic Zeeman coupling qs is infinite, 
the magnitude of F{t) oscillates as the time varies. Prac- 
tically, in most cases we examine below, m^fl turns out 
to be much less than unity because both the curvature 
V^j: and the Zeeman coupling qs are smaller than the in- 
verse of effective mass , or the spin interaction energy 
cl; thus oscillations of the imaginary part of F{t) can be 
neglected. Oscillations along the direction of there- 
fore always lead to temporal oscillations in population of 
rubidium atoms in either |2,0 > or |2,±2 > states that 
can be observed in experiments. 

We first examine the cases when the quadratic Zee- 
man coupling is absent i.e. qB — 0. The frequency in 
this limit is a direct measure of quantum fluctuations 

and Qa — = 0,V) — 2^V^^\cl\- In the absence 

of lattice potentials, or in traps we find that oscillation 
frequencies are about 10~'^Hz. At finite temperatures, 
thermal fluctuations do enhance oscillation frequencies 
by a factor of two to four; however, this is far from suf- 
ficient for the experimental study of QFCSD within the 
life time of these isotopes 10]. 

In optical lattices, we find that the enhancement of 
spin-dependent interactions 6l,cl, and especially the 
rapid increasing of band mass msN can result in oscil- 
lation frequencies of order of a few Hz, which are about 
three to four orders of magnitude higher than those in 
traps. The scaling behavior of VIq is closely related to 
that of Bqf. Taking into account the expression for ef- 
fective mass TO^, we find 
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FIG. 7: (Color online) Frequencies / = n/(27r) for popula- 
tion oscillations as a function of quadratic Zeeman coupling 
gfl(in units of nfc) for various potential depth V in 2D opti- 
cal lattices with M atoms per lattice site, (a) for oscillations 
along the X^-direction around state |2, >; dynamics are un- 
stable above a threshold qci{V). The V-dependence of gd 
(in units of nk) is shown in (b). (c) is for oscillations around 
states with ^ = ^m(< 7r/2) at which global potential minima 
shown in FiglSja) are found, (see also Eq. (|15p ) 



CL 



M\cl\ 

tL 



(17) 



In Fig.m we show the plot of /o = rio/27r, the oscillation 
frequency versus V in the absence of Zeeman coupling 
qs- 

Experimental studies of QFCSD can also be carried 
out by investigating dynamics of rubidium atoms in the 
presence of finite quadratic Zeeman coupling qB- The 
submanifold dynamics now are driven by both quadratic 
Zeeman coupling and quantum fluctuations. There are 
three main modifications to the effective Hamiltonian 
i/g/ / ■ Firstly, according to Eq. ([13]) , effective masses now 
depend on the quadratic Zeeman coupling gg. Secondly, 
following discussions in the appendix^ among four spin 
collective modes [v = x,y,z,^) and one phase mode 
(ly ~ p), two of spin modes, x- and y-mode, can be gapped 
and one, z-mode, always remains gapless. Details of spec- 
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FIG. 8: (Color online) (a), (b) Frequencies / = Q,/{2-k) for 
oscillations along the X^- direction around biaxial nematic 
state |2,2 > +|2,— 2 > versus quadratic Zeeman coupling 
qB (in units of nk) for various optical potential depth V in 
2D optical lattices with M atoms per lattice site, (a) is for 
population oscillations close to the threshold qc2 and (b) is 
for away from qc2- Oscillations are unstable below a threshold 
qc2', in (c), qc2 (in units of nk) as a function of V is shown. 

tra of ^-mode depend on the value of ^; excitations are 
unstable when is negative and are gapped when are 
positive. The energy gap of stable ^-mode excitations de- 
pends on ^ and becomes maximal at ^ = 7r/2. Finally, a 
mean field quadratic Zeeman potential that depends on 
and a = x,y,z, is also present and modifies the 
spin dynamics. Oscillation frequencies O as a function 
of quadratic Zeeman coupling qs for various potential 
depth V can be calculated explicitly using Eq. (fTB|) . Plots 
of these results are shown in Fig |7l8l 

In the limit of weak quadratic Zeeman coupling, 
QFCSD can be most conveniently studied around the 
uniaxial nematic state at ^ = 0, or |2,0 > state 'ss']. 
For rubidium atoms, this is the ground state when the 
quadratic Zeeman coupling is absent, and remains to be 
locally stable along the direction of up to a finite 
Qci. When the quadratic Zeeman coupling (/s is much 
smaller than qd, the effective potential Ve// is mainly 
due to the quantum- fluctuation induced one, Vqf. Above 
gel, a dynamical instability occurs and a perturbation 
along the X^-direction around |2,0 > starts to grow ex- 
ponentially. We find that in traps where V = 0, qd is 
about 10~^pfc; in optical lattices when V varies between 
2Eii to lOEji, the value of qd increases from O.lpfc up to 
about O.lnfc. And as qs approaches zero, the oscillation 
frequency Q{qB,V) saturates at the value Qq shown in 
FiglH The gs-dependence of frequencies for oscillations 
around ^ = is shown in Fig[71 near qs — qd, the fre- 
quencies scale as y/qd — Qb- Alternatively, one can also 
study oscillations around states at which global poten- 
tial minima of V^ff are found; similar g^-dependence is 
shown in FiglTj:. 

At relatively high frequencies, convenient oscillations 




qgLnK] 

FIG. 9: (Color online) Regions of dynamical stability in the 
V — qB plane for three particles per lattice site (M = 3.0). 
Oscillations around uniaxial nematic or state |2, > (biaxial 
nematic or state |2, 2 > +|2, —2 >) are locally stable along the 
X^-direction in the yellow shaded (purple patterned) area; dy- 
namical instabilities occur when gs > qci{V) {qB < qc2{V)). 
qB is in units of nk. 

to investigate are the ones along the X^ direction around 
the biaxial nematic state at f = 7r/2 point or (|2,2 > 
+ |2, —2 >)/\/2 state. This state becomes stable when qs 
is larger than qc2 and oscillations are well defined only in 
this limit. When qB is smaller than this critical value, 
the dynamics are mainly driven by quantum fluctuations 
and oscillations are unstable. In the vicinity of qc2, the 
oscillation frequency again scales as y/qs — qc2- When 
is much bigger than qc2 but smaller than cl, the fre- 
quency is proportional to ^/qs] in this limit, the potential 
Vef / is already dominated by the quadratic Zeeman term, 
or Vq, but the effective mass is mainly induced by spin 
dependent interactions. 

When the quadratic Zeeman coupling further increases 
well above the value of cl, relatively fast dynamics are 
now mainly driven by the external coupling and the 
frequency ft approaches Sgs/S, which is equal to the 
quadratic Zeeman splitting between state |2,0 > and 
|2, ±2 >. In this limit, scattering between atoms during a 
short period of 2w/^l becomes negligible; the population 
oscillations due to scattering or interactions are there- 
fore significantly suppressed as here hyperfine spin two 
atoms are effectively noninteracting. Indeed, following 
the expression for F{t) in Eq. p3)l and Eq . . (fTCj) . one 
finds that the amplitude of population oscillations scales 
as 1/qB at large qs limit and F{t) becomes exp(i4(7Bi/3) 
when qs approaches infinite, i.e., a pure phase factor with 
a constant modulus. 

In FiglSl we show that the gs-dependence of oscilla- 
tion frequencies both close to the threshold qc2 and away 
from qc2- As we have mentioned before, QFCSD is rather 
sensitive to the variation in V . qd and gc2 as functions 
of V are shown in Fig. I7I8I 

So far we have studied coherent dynamics driven by 
potentials induced by quantum fluctuations at zero tem- 
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perature. In the following section, we analyze the effect 
of temperatures, or thermal fluctuations. 

IV. EFFECTS OF THERMAL FLUCTUATIONS 

We now turn to the effect of temperatures and focus on 
three-dimensional optical lattices. At a finite tempera- 
ture, collective modes including spin waves are thermally 
excited. The occupation number of ath spin-wave exci- 
tations with momentum Hq and energy Ect q is 

^('^'1) = ^ • (18) 

e kT — 1 

The free energy density (or per lattice site) of a conden- 
sate characterized by £, is 



Here we have neglected the term depending on T/Tbec 
because we are interested in temperatures much less than 
7B£;c-(Such a term is only important at very low temper- 
atures that are of an order of |cl|(|cl|/Is£;c)^/^ (< T*)). 
Therefore, the low temperature enhancement of Bth is 
mainly from the " black-body" radiation of spin wave ex- 
citations. 

When T T^,, spin modes with energy much bigger 
than are thermally occupied. For modes of energy 
T*, we can approximate n(a,q) with a classical result 
n(a, g) = kT/Ea.q. Our scaling analysis of the free en- 
ergy per lattice site shows that in this limit, 



Fi^,T,V)^—Y, 



1 



a,q 



E 



kT\n{l - e 



.(19) 



Just as in the zero temperature case because of the 
^-dependence of spin wave velocities, the free energy 
density also depends on the values of ^ and the main 
contribution to the ^-dependence of free energy density 
is again from fluctuations of wavelength l/{mBNVa), 
or of a characteristic energy that is of an order of 
\cl\ (or 6l). Thermal fluctuations become more impor- 
tant than quantum ones when temperatures are much 
higher than T*. The other temperature effect is from 
the temperature dependence of spin wave velocities Va- 
Becuase of the thermal depletion, spin wave velocities 
Va{a = x,y,z) decrease as temperatures increase and 
become zero when the BEC temperature Tbec is ap- 
proached. Therefore the temperature dependence of free 
energy are determined by two dimensionless quantities, 
T/T* and T/Tbec- In all cases we are going to study 
below, (of order of — cl) is less or much less than 
Tbec- 

We study the asymptotics of free energy in low temper- 
ature (T <^ Tt.) and high temperature (T 3> T*) limits. 
When T <ti T,,, only modes with energy much smaller 
than are thermally occupied. Therefore contributions 
from these thermal excitations can be calculated by sub- 
stituting quasi-particle spectra Ea^q approximately with 
phonon-like spectra VaQ- One then obtains the following 
expression for the free energy per lattice site (up to a 
constant) 



F{tT,V) = Vqfil 
TT^kT 



Tbec 



90 



E 



dLkT 



(20) 



The barrier height, Bth{T,V) = ^ \.T,V) - 

F{t, = 0,T, y), can be calculated accordingly and the 
asymptotics at low temperatures is 



T 



Tbec 



)^)^.(22) 



The barrier height, Bth{T,V) = F{i = f,T,y) - 
F(^ = 0, T, V), can be calculated accordingly and asymp- 
totics at high temperatures are 



Bth{T, V) 



cl'^ kT 



+3/2 



(l-( 



EEC 



(23) 



The above result shows that at relatively high temper- 
atures thermal spin wave fluctuations enhance induced 
potentials and this enhancement is characterized by a 
linear function of fcT/|ci| or kT/T^. However, because 
of the thermal depletion of condensates at finite temper- 
atures, there is an overall suppression given by a fac- 
tor 1 — {T /Tbec)'^^'^ ■ The competition between these 
two effects results in a maxima in the free energy den- 
sity. That is at temperatures larger than T, but much 
smaller than TbeCi the free energy increases linearly as 
a function of temperature until the thermal depletion be- 
comes significant. The maximal enhancement is therefore 
about Tbec/T^. At further higher temperatures, the 
free energy density decreases and near Tbec, its tem- 
perature dependence is mainly determined by a factor 
(1 - (T/TBEcf^^?^'^ as shown in Eq.^. 

Evidently, the magnitude of maximal enhancement is 
very much dependent on the ratio between Tb ec and . 
In optical lattices, the magnitude of BEC transition tem- 
peratures for averagely one or two particles per lattice is 
mainly set by the band width while is determined 
by \cl\. This ratio therefore the enhancement is quite 
sensitive to the optical lattice potential depth V. We 
find that the enhancement in the absence of lattice poten- 
tials is substantial and thermal fluctuations increase the 
induced potential by about a factor of flfty. This is con- 
sistent with a previous calculation [s^. However, when 
the lattice potential depth increases, increases but 
decreases, the thermal enhancement becomes much less 
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FIG. 10: (Color online) Induced free energy barrier height as 
a function of temperature for different lattice potentials. The 
dashed line in (b) is a plot for the high temperature asymp- 
totics derived in Ea. (|22|) . Note that thermal enhancement of 
barrier height is much less significant in deep optical lattices. 



significant. At = lOEp, we find that the thermal ef- 
fect only increases the potential by a few tens of percent 
and the potential can be predominately due to quantum 
fluctuations. 



V. EFFECTS OF SPIN EXCHANGE LOSSES 

For F ~ 2 rubidium atoms, the life time is mainly 
due to spin exchange losses that take plac e at a rel- 
atively short time scale of about 200ms|Tol ITl| . Note 
that three-body recombination usually occurs at a much 
longer time scale and has little effect on the time evo- 
lution studied here. Spin exchange losses put a very 
serve constraint on possible observations of full coher- 
ent oscillations driven by quantum fluctuations at fre- 
quencies well below 5Hz. To overcome this difflculty, we 
suggest to apply a quadratic Zeeman coupling close to 
<Zci {Qc2) and study time evolution of state |2,0 > or 
^(|2,2 > -|-|2,— 2 >). Regions of dynamical stability 
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FIG. 11: (Color online) Estimated time evolution of rubidium 
population at |2, > hyperfine spin state when the quadratic 
Zeeman coupling or Zeeman field varies, (a) and (b) are for 
optical lattices with potential depth V — WEr and (c) is for 
V — 6Er. In (a), we assume the density is uniform and there 
are no spin losses while (b) and (c) are for traps with inho- 
mogeneous density and a finite loss rate 1/r = (200ms)~^. 
In (a),(b),(c), magnetic fields which yield dynamical critical 
quadratic Zeeman coupling gci(~ qc2) are 35mG, 35mG and 
lOmG respectively. 



of these states can be obtained by studying small oscil- 
lations around them as shown in section III. Results are 
summarized in Fig [51 

Dynamical stabilities of |2, > state below the thresh- 
old qci, or instabilities of (|2,2 > -|-|2,-2 >)/\/2 state 
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below threshold qc2 are directly induced by quantum 
fluctuations. Coherent quantum dynamics around state 
|2,0 > ((|2,2 > +|2,-2 >)/V2) below and above dy- 
namical critical coupling qci{qc2) are qualitatively dif- 
ferent and it is therefore plausible to probe these dy- 
namics that are driven mainly by quantum fluctuations 
at relative low frequencies. To explore this possibility, 
we have studied strongly damped population oscillations 
around a condensate |2, > or a uniaxial state, at differ- 
ent quadratic Zeeman couplings when the spin exchange 
loss time is set to be 200ms. For this part of calcula- 
tions we also take into account the density inhomogene- 
ity of rubidium atoms in traps. Main results are shown 
in FigHT] for different magnetic fields; at a given field, 
the corresponding quadratic Zeeman coupling is given as 
QB = 3{iibB/2)^ /6.8GHz. 

Our major findings are three-folded and outlined be- 
low. Firstly, for a condensate with small deviations from 
|2,0 > state, because of fast spin exchange losses com- 
plete coherent oscillations around |2,0 > are hard to ob- 
serve. Secondly, dynamics can be dramatically modified 
by a finite quadratic Zeeman coupling. When the cou- 
pling is much less than qd (defined fo the density at 
the center of a trap), dynamics become slower when the 
quadratic Zeeman coupling increases. Beyond qd, dy- 
namics become faster while the quadratic Zeeman cou- 
pling increases. These behaviors are consistent with re- 
sults for the case of uniform density shown in the previous 
section; there, the oscillation frequency becomes zero at 
a dynamical critical coupling. 

Thirdly (and perhaps most importantly), the popula- 
tion of atoms at state |2,0 > always grows initially at 
short tome scales when the quadratic Zeeman coupling 
is zero or smaller than qd, as a direct consequence of dy- 
namical stabilities of a |2,0 >-condensate. By contrast, 
beyond qd the population of atoms at state |2, > al- 
ways decreases initially at short time scales as a result 
of dynamical instabilities in this limit. These two dis- 
tinct short-time behaviors of population at state |2,0 > 
are signatures of a transition from quantum-fiuctuation 
driven dynamics to a mainly quadratic Zeeman coupling 
driven dynamics. In the presence of strong spin losses, 
observing these distinct short-time asymptotics is an ef- 
fective way to probe many-body quantum fiuctuations in 
condensates we have considered. At very large quadratic 
Zeeman coupling, we have found expected rapid oscil- 
lations corresponding to mean field coherent dynamics 
and we don't show those results here. Finally, we also 
estimate the time-dependence of population oscillations 
at finite temperatures using thermal-fiuctuation induced 
potentials discussed in a previous section. Thermal fluc- 
tuations usually speed up the coherent dynamics around 
state |2, > when the quadratic Zeeman coupling is suf- 
ficiently weak (smaller than qd)] on the other hand, 
they can also slow down coherent dynamics around state 
|2,2 > -|-|2,— 2 > when a quadratic Zeeman couphng 
larger than (7^2 is present. These results are shown in 
FigHH In obtaining these results, we have extrapolated 
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FIG. 12: (Color online) Estimated time evolution of popu- 
lation at |2,0 > state at different temperatures, (a) is for a 
magnetic field of 13mG that is below the dynamical critical 
field corresponding to qd and (b) is for 38mG field which is 
above qd. The optical lattice potential depth is 1/ = SEu, 
Tbec ~ 10.9nfc and the spin exchange life time is again set 
to be 200ms. Density inhomogeneity in traps has also been 
taken into account to obtain these plots. 

the linear dynamical analysis into a regime where os- 
cillation amplitude is substantial. So strictly speaking, 
results for 35mG, AAmG, 65mG magnetic fields shown 
in Fig lllf b) are qualitative. However, short-time asymp- 
totics that are the focuses of our discussions here are 
accurate. 



VI. CONCLUSION 

In conclusion, we have demonstrated that unlike the 
usual mean-field- interaction driven coherent dynamics, 
QFCSD is a novel class of coherent dynamics fully driven 
by quantum fiuctuations. These dynamics are conve- 
niently tunable in optical lattices where oscillation fre- 
quencies can be varied by three to four order of mag- 
nitude. Frequencies of oscillations driven by quantum 
fluctuations have both distinct lattice-potential depen- 
dence and quadratic Zeeman-coupling dependence that 
can be studied experimentally. For rubidium atoms, 
QFCSD can be directly probed either at frequencies of a 
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few Hz or when the quadratic Zeeman couphng is about 
lO^^nk. One of potential applications of QFCSD per- 
haps is the possibility of performing precise measure- 
ments of strongly correlated quantum fluctuations and 
critical exponents near quantum critical points. Cur- 
rent studies of critical correlations are based on analyzing 
statistics of interference fringes [36l. [37|. Given the great 
control of coherent dynamics recently demonstrated for 
cold atoms [sl], we would like to believe that quantum- 
fluctuation controlled dynamics might be an alternative 
and promising path towards probing critical correlations. 
Finally, we have also investigated the effects of finite tem- 
peratures, spin exchange losses and studied time evolu- 
tion of condensates mainly driven by quantum fluctua- 
tions. 
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APPENDIX A: COLLECTIVE COORDINATES 

To study QFCSD, we expand •0^ about a reference con- 
densate X, 



in terms of ^^(0), 9^(0) for q = modes. 



X'x — 



01 -0p ^ _ e] + 6t 



2WN 



2VN 



01 + Ox 



X,i 



-i^/N iXyy - Xzz)' 
01+ Oy 



WN{xz. - Xxx) 

ol + o. 



WN{xxx - Xyy) 



(A2) 



where = NtM is the total number of atoms in a lat- 
tice with Nt sites. By comparing Eq. (|Aip and the ex- 
pression for X in terms of five spin coordinates X^, (see 
discussions before Eq.®), one identifies the semiclassi- 
cal collective coordinates A^ as the expectation value of 
harmonic oscillator operators A^. In the following, we 
do not distinguish between X^ and X^. The conjugate 
momentum operators P^, can be introduced accordingly 
so that the usual commutating relations are obeyed, i.e 

Quantum-fluctuation-induced potentials Five 
modes of zero point motions (microscopic) around xi^,) 
can also be labelled by the same set of indices. For collec- 
tive excitations of mode-i/ with lattice momentum (\{^ 
0), creation operators are 6'J(q). These operators obey 
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the bosonic commutation relations. 

Note that only the properties of three rotation modes and 
deformation mode (^-mode) depend on the value of ^. 
The collective mode dispersion around ^ = or ^ = 7r/2 
has a very simple form. 



+ -^^ij;^(^t(0) + ^e''i--'=^t(q))(.Al) 

Here Nt is the number of lattice sites and M is the num- 
ber of atoms per site. The superscripts (or subscripts) 
V = x,y,z,^,p specify condensate motion along the five 
orthogonal directions of the submanifold. 9u{<i) (^i(q)) 
are annihilation (creation) operators of the collective ex- 
citations of mode v with lattice crystal momentum q. 
One can further include bilinear terms in the expansion 
(not shown here). The motion along these five orthog- 
onal directions is studied by introducing five mutually 
orthogonal matrices L'^ , i.e. Tr(i^L'^) — 25^^. Ex- 
plicit forms of these orthogonal matrices were introduced 
previously [13. 

We are restricting ourselves to the dynamics of a con- 
densate in a linear regime. By expanding the original 
Hamiltonian in Eq.((T]) using the decomposition intro- 
duced in Eq. ljAip . we obtain an effective Hamiltonian 
that is bilinear in terms of 9'^ and 9. It is convenient to 
introduce five collective Hermitian operators X^ defined 



i;,,q(e,A5,(7B) = ^ tuA'^mBNvli^ + Aj) + e,,q]. (A3) 

And in the expression for £'iy,q, ei/,q = 4ti^^(l — 
co&qadL,) + Ku<lB IS the energy of an atom with crys- 
tal quasi- momentum q = (qx,qy,qz)', cIl is the lat- 
tice constant, tubn — l/(4iLdi) is the effective band 
mass. Kjy, ly — x,y,z,£, are coefficients introduced in 
Eq. p2l) . (|13p . is the sound velocity of the j/-mode 
in the small-|q| limit, f^(C) = M {ibLGaa — cl) /msN, 
w| = M\cL\/mBN- Gapii) is a function of ^ introduced 
in Eq.mkl). 

When the quadratic Zeeman coupling is zero, the en- 
ergy dispersion was also derived in a previous work of 
ours[33|. The dispersion of the ^-mode in this case 
is independent of ^ or A^. And among four gap- 
less spin modes, only the x-, y- and z- spin rota- 
tional modes contribute to the ^-dependence of the 
quantum fiuctuation- induced potential Vqf. Substitut- 
ing Eq. (|A3p into Eq. fT^ . one obtains the quantum- 
fiuctuation-induced potential Vqf. In 3d optical lat- 
tices, we find that the ^-dependent energy potential is 
proportional to "'^a/C'^i^i)' '^l is the lattice con- 
stant and is the hopping integral. In 2d lattices. 
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the ^-dependent energy potential Vqf{^) is proportional 
to Iniidi/wc,. These calculations lead to 

the scaling behaviors discussed in Eg. (fT4|) . When the 
quadratic Zeeman coupling is present, we find that for 
f = 0, Kx,y,z,^ are not positive defined and collective spin 
modes can be unstable. To calculate induced-potential 
Vqf, we include stable modes (q 0) which contribute 
to the renormalization of adiabatic condensate dynamics. 
At ^ = 7r/2, all collective modes are stable. 

Condensate depletion To study the condensate de- 
pletion in the optical lattice, we start from the effective 
Hamiltonian. Five collective modes on top of a nematic 
state are described by the following: 



(A4) 



Here O]^ ^ and 9^^ ^ are bosonic creation and annihilation 
operators for the v-i\i mode. These operators can be 
obtained by Bogoliubov transformations: 



y.q y.q y.q ' y.q u^ — q' 



The coefficients Uy^q and Vy^q are given: 




1/2 



1/2 



(A5) 
(A6) 



(A7) 
(A8) 



Condensate deletions from the v-ih. mode with crystal 
momentum q are 

{Olj.^q) = Kq + Vlq) {O^qK^q) + (A9) 

Here {dl^qO^ g) is the occupation number of the vth mode 
and it is given by the Bose-Einstein statistics: 



Putting all these together, we reach the final equation to 
determine the condensate density: 



M - Mo{T) 




(All) 



M is the atom number density or number of atoms per 
lattice site and Mq is the condensed number density. 
The first term in the summation gives quantum depletion 
number due to the two body scattering and the second 
term gives the thermal depletion number. We solve the 
above equation numerically and the results are shown in 
Figffl 

At zero temperature, quantum depletion is given by 
the term ^ which captures the effect of interactions or 
two body scattering processes. Our estimate shows that 



Mo ^ V — 



1 (MaL 



M \ tLdL 
3/2 



(A12) 



We have noticed that major contributions are from the 
phase fluctuations since 3> |ci|. In optical lat- 
tices when the potential depth V increases, the bandwith 
decreases exponentially as a function of V and az, in- 
creases; as a result, quantum depletion grows as {y^Y^"^- 

At finite temperatures, Mq decreases as temperature 
T increases and becomes zero at the transition temper- 
ature Tbec- In the weakly interacting limit (t^ >> 
olj^LiIclI), the condensate fraction is approximately 
equal to 



Mo(T) 
Mo(r = 0) 



= 1 - 



T 



BEC 



(A13) 



exp {E^,q/kT) - 1 ' 



(AlO) 



2 /3 . 

and Tbec ^ 5Mq in optical lattices. 
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